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We compute the contribution of blazars (Flat Spectrum Radio Quasars and BL-Lacs) to the 
Fermi-LAT Extragalactic 7-ray Background, assuming a 7-ray luminosity function proportional to 
the luminosity function of radio galaxies 9] and the Spectral Energy Distribution of blazars given by 
the so-called blazar sequence We find that blazars account for the ~45 % of the total Fermi-LAT 
EGB. However, at low (<10 GeV) and high (>50 GeV) energies the blazar contribution falls short 
to account for the observed background light. At E<10 GeV, the Extragalactic 7-ray Background 
can be fully explained by considering the additional contribution of star-forming galaxies, whereas 
at very high energies, where the absorption due to Extragalactic Background Light dominates, the 
observed background calls for an additional unknown component. 



INTRODUCTION 



The Extragalactic 7-ray background (hereafter 
EGB) represents a fascinating challenge since its first 
detection by SAS-2 satellite above 30 MeV [J. The 
7-ray telescope EGRET during nine years of opera- 
tions, measure an isotropic 7-ray emission in the 30 
MeV-30 GeV band @. The satellite Fermi-LAT after 
one years of observations improves the EGRET mea- 
surement confirming the existance of an EGB up to 
100 GeV [|. 

As already found by EGRET and confirmed by 
Fermi-LAT 0], blazars are the most common ob- 
jects in the extragalactic 7-ray sky. Blazars, radio 
loud AGNs with the jet pointing towards the Earth, 
are classified as Flat Spectrum Radio Quasars (FS- 
RQs) and BL-Lacs. After three years of observations, 
Fermi-LAT detected approximately 1000 blazars with 
about the same number of FSRQs and BL-Lacs 
It results to be natural computing their contribution 
to the EGB (e.g. The main ingredients in 

the EGB calculation are the 7-ray luminosity func- 
tion (LF) of blazars and their spectral energy distri- 
bution (SED). In our model adopt a LF of FRI and 
FRII radio galaxies Q and the blazar SEDs using the 
the so-called blazar sequence [lol IllT ] . We compute the 
blazar contribution to the Fermi-LAT EGB fitting the 
differential Fermi-LAT logN-logS, having as a sole fit- 
ting parameter the fraction of radio galaxies beamed 
towards the Earth, hence appearing as blazars. Bca- 
cause the blazar component alone is not able to fully 
account for the EGB at low energies (<10 GeV), we 
consider the 7-ray emission of star-forming galaxies. 
This allows us to fully explain the EGB but the last 
point at approximately 80 GeV where the EBL ab- 
sorption strongly dominates. Notice that here we refer 



to EGB as the sum of the resolved source component 
and the unresolved and/or truly diffuse emission. 



II. BLAZAR CONTRIBUTION 



The blazar contribution (in photons s 1 cm 2 sr 
MeV ) to the EGB at the observed energy Eq is 
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where <i<I> 7 (L 7 , z) / d\ogL 1 is the 7-ray LF and L 7 is 
vL v (in erg/s) at 100 MeV, dn(L 7 , z)/dE is the un- 
absorbed photon flux per unit energy E = Eq(1 + z) 
measured on Earth of a blazar with luminosity L 7 at 
redshift z, and r 77 (i?o, z) is the optical depth for 7 — 7 
absorption. We adopt the Extragalactic Background 
Light (EBL) model by [HI]. In the above equation 
dV/dz is the comoving cosmological volume. We set 
logi™ = 43.5 and logi™* = 50. 

The number of sources N(> -Fph) per steradian with 
photon flux greater than F p h is 



N(> F ph ) = i- 



dz — — 

dz 



log L" 



diogL 



(L~f,z) 



d log L 1 
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The 7-ray LF of blazars is presently uncertain (see e.g. 
(HI), so that one has to rely on the LFs computed in 
other bands, e.g., X-rays (i. 14|, or radio 0. We adopt 
here the radio LF at 151 MHz of FRI and FRII [g, 
assumed to be the parent populations of blazars: 



$ 7 (L 7 ,z) <S> R (L R ,z) 



d log L 7 



d log L R 



(3) 
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FIG. 1: The black line represents the fit to the blazar 
logN- logS measured by Fermi [l^| as a sum of BL-Lac 
component (blue line) and FSRQ component (red line). 
Black points are the number counts of all the blazars, red 
points FSRQs and blue points BL-Lacs 



FIG. 2: The black points represent the total EGB (re- 
solved and unresolved sources and all the rest), blue line 
is the FSRQ component, the green line is the BL-Lac com- 
ponent. The red line is the sum of the two components. 



where Lr, is vh v at 151 MHz, and the constant k is 
the fraction of blazars over all radio galaxies. In order 
to convert radio into 7-ray luminosity, we rely on the 
blazar spectral energy distribution (SED). As pointed 
out in [ll(, a relation between the radio and/or bolo- 
metric luminosity and the energy of the two peaks ex- 
ists (the so-called blazar sequence, see also [HI). We 
use the SEDs computed by d, [l(| • 



A. Results 

We use the model discussed in the previous section 
to compute the total contribution of blazars (FSRQs 
and BL-Lacs) to the Fermi-LAT EGB. To this end we 
fit blazar logN-logS distribution by letting k to vary. 
We find k = (3.93 ±0.01) x 10~ 4 . As shown in Fig. [TJ 
once fitted the total differential logN-logS, our model 
predict the correct number of BL-Lacs and FSRQs 
without any change in the radio LF. This fact sup- 
port the correctness of the initial hypothesis of radio 
galaxies as parental population of blazars. The num- 
ber ratio of blazars to radio galaxies k can be thought 
as a measure of the beaming factor of the relativistic 
jet, which in turn is related to the bulk Lorentz factor 
T. From k ~ 1/2T 2 we derive T ~ 35. 
Fig. [5] shows the corrisponding contribution of FS- 
RQs and BL-Lac to the total Fermi-LAT EGB. 
The background intensity is found to be Ifsrq = 



4.22 x 10~ 6 ph s^cm^sr -1 and I B L-Lac = 2.43 x 
10 _6 ph s~ 1 cm~ 2 sr _1 for FSRQs and BL-Lacs respec- 
tively. The total EGB intensity is therefore I = 
6.65 x 10~ 6 ph s _1 cm _2 sr _1 , corrisponding to 45% 
of the one measured by Fermi-LAT . 
From the slope of the FSRQ and BL-Lac component 
in Fig. [5] we can see that the main photon index of 
FSRQs and BL-Lacs resulting from our model are in 
agreement with the 2LAC Q. 

We note that blazars fall short to explain the mea- 
sured EGB at E<10 GeV and at E>50 GeV. At low 
energies, the discrepancy can be fully accounted by 
star-forming galaxies modeled following the recipes 
by[7|, so that only the last point of the Fermi-LAT 
EGB measurement is not reproduced by our model. 



III. DISCUSSION 

We have computed the overall contribution of 
blazars to the Fermi-LAT EGB. 

Our model relies on two assumptions: the radio LF 
and the blazar SED. In the following we show the 
difference with the most reliable works on the contri- 
bution of blazars to the Fermi-LAT EGB. 

• We use the radio LF . Fitting the differential 
logN-logS we obtain just the overall normaliza- 
tion k, without any change on the bright and 
faint end of LF. 
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FIG. 3: The star-forming galaxy component (green line), 
the total blazar component (red line) and the sum (blue 
line). In black are the Fermi-LAT EGB points. 



In previous works [6] 14] it is assumed a LF in X 
band with three free parameters: the total nor- 
malization, the amount of bolometric radiation 
emitted in X-ray and the faint end of the X-ray 
LF. Differently in Q the radio LF computed by 
((l6j) is used, changing the faint end to obtain 
the contribution of FSRQs. 

The second main point of our model is repre- 



sented by the blazar SED. We assume that the 
spectra of blazar is fully taken into account by 
the blazar sequence On the contrary in 

@, [1] it is assumed the blazar spectra as a sim- 
ple and broken power law, respectively. In these 
works they assume a spectral index distribu- 
tion peaked on the mean spectral index resulting 
from observations. The underlying assumption 
is that the unresolved blazar component has the 
same index distribution of the resolved compo- 
nent. Using the blazar sequence such assump- 
tion is not necessary because the blazar SEDs 
are fully determined. 

The best fit value of the relative number of blazars 
with respect to radio galaxies can be translated into a 
bulk Lorentz factor of the relativistic jet T ~ 35, larger 
than the average value V ~ 15 estimated in [la ]. The 
two values could be reconciled if blazars commonly 
show secular 7-ray large variability which modulates 
the 1-year average flux, as recently proposed in [l7| . 
To be consistent with the Fermi-LAT points at lower 
energy, we add the star-forming component to our 
blazar model. Fitting the Fermi-LAT EBG with this 
two component model (Fig. [3]) , we constrains the so- 
called " star formation efficiency of molecular hydro- 
gen" £ = 7.0 x 10~ 10 yr _1 , which we found well within 
existing, much looser observational constraints [l9j . 
The point at 80 GeV remains unaccounted by our 
models calling for an additional (local or low-z) com- 
ponent. An intriguing possibility is represented by 
annihilation of a dark matter particle. We found, as 
an example, that a particle of mass ~ 0.5 TeV and 
cross section (av) ~ 5 x 10 -26 cm 3 s _1 can easly ac- 
comodate the last data point. 
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